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Averaging in general relativity is a complicated operation, due to the general covariance of the 
theory and the non-linearity of Einstein's equations. The latter of these ensures that smoothing 
spacetime over cosmological scales does not yield the same result as solving Einstein's equations 
with a smooth matter distribution, and that the smooth models we fit to observations need not be 
simply related to the actual geometry of spacetime. One specific consequence of this is a decoupling 
of the geometrical spatial curvature term in the metric from the dynamical spatial curvature in 
the Friedmann equation. Ifere we investigate the consequences of this decoupling by fitting to a 
combination of HST, CMB, SNIa and BAO data sets. We find that only the geometrical spatial 
curvature is tightly constrained, and that our ability to constrain dark energy dynamics will be 
severely impaired until we gain a thorough understanding of the averaging problem in cosmology. 



I. INTRODUCTION 

The standard model of cosmology relies on the assump- 
tion that the Universe is well described at all points in 
space and time by a single linearly perturbed Friedmann- 
Lemaitre-Robertson- Walker (FLRW) geometry (except 
in the vicinity of black holes and neutron stars), and that 
this geometry obeys Einstein's equations. However, in an 
inhomogeneous universe this is almost certainly not true. 
Even if the gravitational interaction behaves exactly as 
Einstein predicted on small scales, this will not be true 
for large-scale averages. A critical problem in cosmology 
is therefore determining the form of deviations from Ein- 
stein's equations when considering geometry averaged on 
large scales, and what the effects of these deviations will 
be on observations (for a review, see [T]). This has been 
the subject of some controversy, with opinions ranging 
from the suggestion that the effects of averaging could 
completely explain the recently observed accelerating ex- 
pansion of the Universe without the need for any dark 
energy [2HH], to the claim that it is completely negligible 
[T0H2I] • Others suggest that while the effects of averag- 
ing may not be responsible for the apparent acceleration, 
they may be important for precision cosmology ^22n52] . 

In this paper we take the solutions to the field equa- 
tions derived using an exact and fully covariant averaging 
procedure, and compare them to observations. These so- 
lutions have decoupled spatial curvature parameters in 
the metric and the Friedmann equation, and reduce to 
the FLRW solutions of Einstein's equations when these 
parameters are equal. We find that that the constraints 
available on the spatial curvature parameter appearing 
in the Friedmann equation are considerably weaker than 
those available on the spatial curvature parameter ap- 
pearing of the macroscopic metric. In particular, the con- 
straints from the Cosmic Microwave Background (CMB) 
are considerably weakened, and no longer signal a flat 
universe. This allows for the possibility of averaging hav- 



ing non-negligible dynamical consequences. We also find 
that some data sets prefer models in which the two cur- 
vature parameters are not equal. 

II. SPACETIME AVERAGING AND 
MACROSCOPIC GRAVITY 

There are a number of averaging procedures that have 
been introduced in order to study the large-scale evo- 
lution of inhomogeneous spacetimes. An exact and co- 
variant approach that allows tensor quantities to be av- 
eraged, as well as scalars, was provided by Zalaletdinov 
[33] . Here the geometric objects that exist on the space- 
time manifold are averaged, and the field equations that 
these quantities satisfy are constructed. This is achieved 
using bilocal averaging operators over closed regions of 
spacetime, E, that contain the supporting points x (see 
[33 ) . The result of averaging X is then denoted by {X) . 

Using this definition we can now consider the average 
of various geometric objects. Following Zalaletdinov, we 
denote the average of the connection as (F'^^^^), and de- 
fine a new macroscopic Riemann tensor 

+ (F^„)(F%^)-(F%)(r%J. (1) 

Crucially, M^^^^ ^ {R^u^p)^ ^here (i?''^p^) is the av- 
erage of the microscopic Riemann tensor. From these 
quantities one can then construct the macroscopic field 
equations |33j : 

(/0M,^-i<5%(O^V (2) 

where (^1^) is the averaged energy-momentum tensor, 
{g^^) is the averaged metric, and Z°'^^p = 2Z"^j^''^^j is 



2 



a correlation tensor defined by 
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where underlined indices are not included in anti- 
symmetrization. This quantity must obey the differ- 
ential constraint -^"^[^^jyo-.A] ~ 0, where the covariant 
derivative here is defined using the averaged connection 
(r'^^p). The macroscopic field equations ([2| replace Ein- 
stein's equations on large scales. 



A. Macroscopic FLRW Solutions 

The solutions to the macroscopic field equations 
with an FLRW ansatz for the macroscopic metric, have 
recently been studied in 24, 34 37J, where it was found 
that the extra terms involving the correlation tensor 
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take the same form in the macroscopic field 



equations that a spatial curvature curvature term takes 
in Einstein's equations. In fact, for a spatially flat macro- 
scopic metric, and with spatial correlations only, it can 
be shown that extra terms in Eq. ^ can only take the 
form of spatial curvature. Any other form would be in- 
compatible with either the conservation equations, their 
integrability conditions, or the algebraic constraints that 
^"is[-y'^^a] ni'ist satisfy for consistency of the averaging 
scheme. Thus, the averaged Einstein field equations for 
a spatially flat, homogeneous, and isotropic macroscopic 
spacetime geometry take the form of the Friedmann equa- 
tions of general relativity for a non-flat FLRW geometry. 
That is, the spatial curvature of the macroscopic space- 
time is decoupled from the spatial curvature that appears 
in the macroscopic Friedmann equation. This is an im- 
portant difference from the standard approach to cosmol- 
ogy, where it is assumed that Einstein's field equations 
hold whatever the smoothing scale, and that the spatial 
curvature in the Friedmann equation is therefore identi- 
cal to the spatial curvature of the macroscopic spacetime. 

Let us now write the line-element of the macroscopic 
geometry as 



ds'^ = {g^,^) dx^'dx" = -dt^ + a^{t) 



dr^ 



1 — k^r"^ 
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where the geometrical curvature, kg, is, in general, a func- 
tion of the scale of S. The scale factor a{t) is that of 
the macroscopic spacetime. On scales larger than S the 
macroscopic field equations, ([2]), then become 
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where the 'dynamical curvature', fed, is again a function 
of scale, and we have included in fed here contributions 
from both the spatial curvature in the metric, kg, and 
the terms in Eq. ^ that involve the correlation ten- 



vanishes do we recover the usual result kg = kd, and 
even in this case we can still have a scale dependence. 
More generally, in spacetimes that are inhomogeneous 
on small scales, we do not expect these two 'spatial cur- 
vature' terms to be equal. Defining flk,, — —kg/a^H^ 
and rjfej = —kd/a^H^ the macroscopic Friedmann equa- 
tion then becomes 1 = f2„i + + ^^A, where and 
r^A are the usual expressions for the fraction of the en- 
ergy content of the Universe in matter and the cosmo- 
logical constant, respectively. Note that 17 a is a function 
of smoothing scale whereas A is not. The quantity 0,^ 
does not have to satisfy a constraint of this kind, as it is 
now decoupled from the Friedmann equation. Thus we 
arrive at a parametrized phenomenological model within 
which we can analyze data in order to study the potential 
observational effects of averaging. 



B. Observables in the Macroscopic Universe 

We will now consider distance-redshift relations in the 
FLRW solutions of macroscopic gravity. These provide 
the basis for many key observational tests of the cosmo- 
logical background. 

First we need to know the trajectories of photons in the 
macroscopic geometry. We will take these to be null tra- 
jectories with respect to the macroscopic metric that has 
been constructed to approximate the distance between 
two points in spacetime separated by scales above that 
of S. We consider this to be a reasonable assumption 
for the average of a large number of photon trajectories, 
but note that it will not be true for each individual null 
curve of the microscopic spacetime. If this assumption is 
wrong then it could lead to profound differences with the 
results of applying Einstein's equations directly to non- 
local averaged quantities [351 • Our approach should 
therefore be considered a conservative one. 

Let us now derive the luminosity distance-redshift re- 
lation in the macroscopic geometry. Integrating a null 
trajectory in the geometry (HI), assuming il^^ and Vlk^ 
are constant, and using the solutions to the macroscopic 
Friedmann equation gives 
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where fkgix) = sinh(a;), x or sin(a:) when kg < 0, kg = 
or kg > 0, respectively. This expression reduces to the 
usual one when flf, — flk^. 



III. DATA ANALYSIS 

We can now compare our averaged models to various 
cosmological probes. 

Hubble Rate. A vital tool for constraining dark en- 
ergy is the local Hubble rate. We use Hubble rate 
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data from HST measurements [H], which give Hq = 
74.2±3.6kms~-'^Mpc~^, and we assume Gaussian errors. 

Cosmic Microwave Background. The CMB is well 
known to tightly constrain spatial curvature in the stan- 
dard cosmology. When = fife^ constraints from the 
WMAP 7 year data release [42], combined with HST 
[5T] and Baryon Acoustic Oscillation (BAO) data [33], 
gives -0.0133 < flk < 0.0084 (95% CL). This constraint 
arises principally from the tight bounds on the area dis- 
tance to the surface of last scattering, which must be 
(1 + z^)dA{z^) = 14150 ± 150 Mpc. We fit three param- 
eters which are predicted by the model: The decoupling 
epoch (z*), the acoustic scale (Ia) and the shift param- 
eter (R), which are sufficient to capture the constraints 
from the CMB . To obtain likelihoods, we also use the 
inverse covariance matrix for the WMAP distance priors 
for these parameters is given in Table 10 from |42j . 

Supernovae. Another key probe of the large-scale ex- 
pansion of the Universe is the observation of type la 
supernovae (SNIa). These events are considered to be 
"standardizable candles", in that their absolute magni- 
tude can be approximated when 'stretch' and 'color' pa- 
rameters have been extracted from fits to light-curve tem- 
plates. They then allow the expansion history of the Uni- 
verse to be mapped, and are widely considered to be one 
of the most compelling sources of evidence for the exis- 
tence of dark energy. We therefore consider them here in 
the context of the FLRW solutions to macroscopic grav- 
ity, in order to determine the consequences of allowing 
rikg ^ ri/cj. The supernova data used in obtaining these 
constraints are the Union2 data set [44 , and the SDSS 
data set [Hj. 

Baryon Acoustic Oscillations. Observations of BAOs 
provide a direct measurement of the Hubble rate at non- 
zero redshifts, and are therefore a powerful tool for con- 
straining dark energy. However, the interpretation of 
BAOs relies on assumptions about the evolution of struc- 
ture in the Universe that may not be valid if averaging 
is important. We therefore choose to use BAO data only 
sparingly. We use the fraction of the comoving sound 
horizon to volume distance for the two points, at red- 
shifts z = 0.2 and 0.35. The inverse covariance matrix is 
given by equation (5) from j43| . 



A. Parameter Constraints 

We use the Monte-Carlo Markov Chain (MCMC) 
method to obtain the marginalized errors on the model 
parameters from the likelihood function, using the pub- 
Ucity available package COSMOMC [iO]. 

First consider CMB-|-iJo constraints, as these lead to 
extremely tight constraints on spatial curvature in the 
standard model. It can be seen from Fig. [T] that when 
7^ the CMB -I- Hq no longer constrains spatial 
curvature significantly, due principally to a degeneracy 
between the effects of flk^ and fi^^ in Eq. (l6|. In fact, 
even with A = there exist values of fcg and that sat- 
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FIG. 1: Constraints on fi^^ and Qk^ from the CMB for 
Ho — 50, 70 and 90 km s~^ Mpc~^, with the other parameters 
marginalised over. Shaded areas are 68% and 95% confidence 
regions. 

isfy the observations. These results significantly weaken 
a key part of the evidence for both a spatially fiat uni- 
verse and a non-zero value of A. 

In Fig. [2] we show the combined constraints that can 
be imposed on il^^ and Qk^ using data from the HST, 
the WMAP 7 year data of the temperature-temperature 
correlations in the CMB, the Union 2 and SDSS SNIa 
data sets, and the constraints on the 'volume distance' 
from the BAOs. Fig. [3] shows the constraints available 
on and Q\ from the same data sets. 

It is clear that the effect of allowing fi^ and il^^ to 
be independent has considerable consequences for these 
probability distributions. The marginalized posterior 
values of each parameter in the various cases are given in 
the table. These values are considerably different to the 
case where flk^ = Vlk^ , which are shown in the same table 
for comparison. It can be seen that the additional free- 
dom gained by allowing 7^ ^k^ is considerable, with 
constraints on J^a and the two fi/c's being significantly 
weaker than in the standard approach. The combination 
of all of these observables, however, still appears to pro- 
vide strong evidence for the existence of dark energy, and 
is still consistent with a spatially fiat universe. Neverthe- 
less, it is striking that constraints on flk^ are an order of 
magnitude tighter than those on fifc^. 



B. Constraints on Dark Energy Dynamics 

Decoupling the geometrical and dynamical spatial cur- 
vature parameters does not appear to allow enough free- 
dom to entirely account for the dark energy component 
(which is now signaled almost entirely by supernova ob- 
servations, rather than by both CMB-|-iJo and SNIa ob- 
servations, as in the standard model). Nonetheless, con- 
sidering averaged spacetimes does considerably weaken 
our ability to constrain dark energy in a meaningful 
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FIG. 2: Constraints on the two curvature parameters from different data sets. In the left-hand panel we show the constraints 
from the CMB (gray) and SNIa (blue for Union2, and hollow and dashed for SDSS) separately, as well as the combined 
constraints including HST. In the right-hand panel we show the combined constraints, with the smaller lightly shaded regions 
in the foreground now including the BAO. The constraints from the CMB and SNIa individually are extremely weak, but 
tighten when combined with HST data. Note that only the geometrical curvature is tightly constrained, and not the dynamical 
one. The SDSS SNIa data can also be seen to inconsistent with Qkg = at greater than 95% confidence. 



Constraints on curvature and A with decoupled curvature parameters, and in the standard model. 
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FIG. 3: Constraints on ilm and Ha for the same data set 
combinations as in Fig. |2] The constraints from the CMB 
and SNIa separately are again much weaker here than in the 
standard case. In particular, A = is consistent with the 
CMB data for any Urn- Union2 SNIa data, however, still 
clearly gives A 7^ 0. 
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FIG. 4: Marginalized constraints on the dark energy param- 
eters, with (red) and without (blue) assuming fcg = fed. Al- 
lowing for effects due to averaging removes our ability to con- 
strain the possible evolution of the dark energy equation of 
state. 
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way. For examples, if we parameterize the dark energy 
equation of state as w{z) = wq + Wa{l — a) = wq + 
Waz/{l + z) then CMB+HST+BAO+SNIa(Union2) give 
only the weak constraint wq = — 1.612^q524, with almost 
no meaningful constraints on Wa at all. This can be com- 
pared with Wo = -1.066tg:i^|7 and Wa = -0.130ti;i35 
when kg = k^, as in Fig.|4j The uncertainty arising from 
the averaging problem is clearly hugely amplified when 
we try and constrain the dynamical properties of dark 
energy. 



C. Scale Dependent Curvature Parameters 

As well as kg and k^ being allowed to be different, 
a second consequence of using an averaged geometry is 
that these parameters could become scale dependent. 
This possibility introduces a considerable amount of ex- 
tra freedom, and hence we consider it separately here 
from the more restricted case we have considered so far 
(where kg and kd are different, but scale independent). 

The cosmological probes we have discussed all make 
observations over different scales. The scale of BAOs at 
the present time is about 150 Mpc. High redshift SNIa 
(out to z ~ 1) cover a range of scales out to several 
Gpc, while the CMB involves making observations on 
the scale of the horizon (~ 14 Gpc). By introducing a 
scale dependence into kg and k^ we can therefore poten- 
tially re-introduce the effects of spatial curvature on the 
scales of SNe and BAOs, while still satisfying the strin- 
gent constraints available on the largest scales from the 
CMB. The question then becomes whether or not it is 
possible for observations to effectively constrain spatial 
curvature on different scales in the Universe. A positive 
detection of conflicting measurements of spatial curva- 
ture on different length scales would be a sure sign of 
non-trivial averaging effects [21 US] ■ 

As we have shown, the CMB may no longer imply spa- 
tial flatness or a non-zero flA at all (see Fig.jl]). If spatial 
curvature on the scale of SNIa observations is allowed to 
be independent from the constraints imposed by other 
observables, then the evidence available for the existence 
of a non-zero is again severely weakened (see Fig. |3]). 
Also, if ilj;^ and fi^^ vary sufficiently rapidly on scales of 
interest for SNIa observations then we can no longer treat 
them as being simple constants when analyzing the data. 
Instead, since we average these observables over the sky, 
what we are doing is effectively averaging the geometry 
out to some redshift. This should be expected to result 
in redshift dependent effective curvature parameters, and 
in this case one could end up with curvature parameters 
that are effectively functions of radial distance, k = k{r). 

This picture is somewhat similar to LTB void models 
of the Universe, where we are at the center of a large 
spherically symmetric inhomogeneity. It is well known 
that such models are able to explain the supernova data 
without evoking the existence of dark energy j46j , but at 
the expense of strongly violating the Copernican Princi- 



ple. In the present interpretation no such violation need 
occur, as all observers would experience a universe with 
k — k{r), with themselves at the center of symmetry. 
This would relieve the key philosophical problem associ- 
ated with these models as an explanation of the data. It 
would also relieve the strong constraints on these models 
that are available from the kinematic Sunyaev-Zeldovich 
effect [47-49j, as every cluster would experience (approx- 
imately) isotropic CMB radiation, just as we do. How- 
ever, such a departure from the usual interpretation of 
LTB models would also undoubtedly require revisiting 
the problem, as the field equations would be modified by 
additional terms due to averaging. 



IV. CONCLUSIONS 

We have presented and constrained models of the large- 
scale Universe that result from averaging the geometry 
of spacetime. These models have decoupled spatial cur- 
vature parameters in the macroscopic line-element (f^fc^) 
and the Friedmann equation (0^^), and provide a qual- 
itative alternative to the standard model of cosmology. 
They can therefore be used to analyze the statistical sig- 
nificance of the standard model in a larger space of mod- 
els that allows for some of the non-trivial consequences 
of averaging. 

Using HST, CMB, BAG and SNIa data it is clear that 
the effect of allowing flk and flk^ to be independent 
has considerable consequences for parameter estimation. 
Analysis of the available data shows that the size of the 
68% and 95% confidence regions of f^A, and ilka are 
all much larger than in the standard model. There are 
even tantalizing hints that the data may favor flk ^ 17^^ 
(the combination of HST, CMB and SDSS SNIa data ex- 
cludes flfc^ = rjfcj at the 95% confidence level). However, 
while the evidence for JIa ^ available from individ- 
ual observables can be considerably reduced, the combi- 
nation of SN, CMB, HST and BAG data still provides 
strong evidence for the existence of dark energy, as long 
as rjfe and fi^^ are scale independent universal constants. 
Relaxing this last assumption makes combining observ- 
ables on different scales a much more complicated prob- 
lem, and it is highly probable that such additional free- 
dom will significantly weaken the evidence for SIa 7^ 0. 
The way that this situation should be modeled, and con- 
strained with data, is still an open problem. 

Finally, we have shown that introducing uncertainty 
due to averaging into our models of the Universe dra- 
matically weakens the constraints that can be imposed 
on the equation of state of dark energy. As a conse- 
quence, it therefore seems necessary to understand and 
incorporate the effects of averaging in general relativity 
into our models if we are to begin attempting precision 
cosmology. 
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